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Abstract. Max cones are the subsets of the nonnegative orthant R^ of the 
n-dimensional real space R" closed under scalar multiplication and componen- 

O. twise maximisation. Their study is motivated by some practical applications 

' which arise in discrete event systems, optimal scheduling and modelling of syn- 

chronization problems in multiprocessor interactive systems. We investigate 
the geometry of max cones, concerning the role of the multiorder principle, 
r^ , the Kleene stars, and the cyclic projectors. 

The multiorder principle is closely related to the set covering conditions in 
max algebra, and gives rise to important analogues of some theorems of convex 
geometry. We show that, in particular, this principle leads to a convenient 
representation of certain nonlinear projectors onto max cones. 

The Kleene stars are fundamental in max algebra since they accumulate 
\l ' weights of optimal paths and yield generators for max-algebraic eigenspaces of 

^ ' matrices. We examine the role of their column spans called Kleene cones, as 

building blocks in the Develin-Sturmfels cellular decomposition. Further we 
Cn i show that the cellular decomposition gives rise to new max-algebraic objects 

Oa ' which we call row and column Kleene stars. We relate these objects to the max- 

(_P ' algebraic pseudoinverses of matrices and to tropical versions of the colourful 

^^^ , Caratheodory theorem. 

^"^ ■ The cyclic projectors are specific nonlinear operators which lead to the 

so-called alternating method for finding a solution to homogeneous two-sided 
^~^ ' systems of max-linear equations. We generalize the alternating method to the 

case of homogeneous multi-sided systems, and wc give a proof, which uses the 

cellular decomposition idea, that the alternating method converges in a finite 

1 J ■ number of iterations to a positive solution of a multi-sided system if a positive 

|^N( ' solution exists. We also present new bounds on the number of iterations of 

the alternating method, expressed in terms of the Hilbert projective distance 
between max cones. 
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1. Introduction 



The nonnegative orthant M" of the n-dimensional real space R" can be viewed 
as an n-dimensional free semimodule over the max-times semiring, which is the set 
of nonnegative numbers IR+ equipped with the operations of 'addition' a (B b :— 
max(a, b) and the ordinary multipUcation a (E) b := a x b. The max-times semiring 
is denoted by Mmax.x = (M+,© = max,® — x). Zero and unity of the semiring 
coincide with the usual and 1. For instance, in this semiring 2 (8) 3 = 6 and 
2 © 3 = 3. Subsemimodules of R" — RJ^j^x x ^^6 the subsets of R" closed under 
the componentwise maximization ©, and the usual multiplication by nonnegative 
scalars. These subsemimodules will be called max cones, due to their obvious 
analogy with convex cones. In a very important special case, max cones can indeed 
be convex cones, but in general they are not convex, i.e., not stable under the usual 
componentwise addition. 

By max algebra we understand linear algebra over the semiring Rmax, x , ex- 
tending the max, x arithmetic to nonnegative matrices and vectors in the usual 
way. For instance, HA— (uij) and B = {bij) are two matrices of appropriate 
sizes, then {A © B)ij = aij © bij, or {A (8) B)ij = 0^, atkbkj- The iterated product 
A^ A® ... ^ A in which the symbol A appears k times will be denoted by A''. We 
assume that A" := /, the unit matrix. The sets like {1, . . . , m} or {1, . . . , n} will 
be denoted by [m] or [n] respectively, and for a set of indices M, the number of 
elements in M will be denoted by \M\. 

The idempotency of addition a(Ba — a and the lack of subtraction are important 
features of max algebra that make it different from the nonnegative linear algebra. 

Max algebra has been known for some time, and we mention here the pioneer- 
ing works of Cuninghame-Green [18L 119] . Yoeli [47) . Vorobyev [45], Carre [12| . 
Gondran and Minoux 27 , K. Zimmermann |48) . and U. Zimmermann [50], among 
many others. Max algebra is often presented in the settings which seem to be differ- 
ent from Rmax,x , namely, over semirings Rmax,+ = (R U {— oo}, © = max, = -I-) 
(max-plus semiring), R,iii,i.+ = (R U {+oo}, © — min, = -f ) (tropical or min-plus 
semiring), or most exotically Rmin.x ~ (R+ U {+oo},© ~ min, = x) (min-times 
semiring). All these semirings are isomorphic to each other and to Rmax,x- Max 
algebra has important practical applications which arise in discrete event systems 
and scheduling problems [2l 1191 123] , and in modelling of synchronization problems 
in multiprocessor interactive systems |10| . 

More generally, max algebra can be seen as a branch of tropical mathematics, 
which is a rapidly developing field with applications in mathematical physics, op- 
timal control, algebraic geometry and other research areas. See |36| for a recent 
survey, and ^34, 35 for recent collections of papers. 

The similarity between max cones and convex cones was understood in the very 
beginning by Vorobyev [45|, who used the name 'extremally convex cones' (instead 
of semimodules or spaces). K. Zimmermann [49j defined extremally convex sets, or 
tropically/max-plus convex sets as it would be called now, and proved a separation 
theorem of a point from a closed convex set. This theorem was generalized and 
more transparent proofs were given by Samborskii and Shpiz [42], Litvinov et al. 
[37) . Cohen et al. [HI [S], and also Devehn and Sturmfels [22], Joswig [32]. We 
note that the separation theorem of a point from a closed max cone, given below 
as Theorem 12. 6[ is essentially the same result. In the ordinary convex geometry, 
separation of a point from a convex set easily leads to the separation of two convex 
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sets from each other. However, analogous statements for max cones arise differently 
and are related to the investigation of certain nonlinear projectors onto max cones, 
and their compositions called cyclic projectors, see Gaubert and Sergeev |26) and 
Theorems 14.21 and 14.31 below. Remarkably, these cyclic projectors also appear in 
the study of two-sided max-linear systems of equations, see Cuninghame-Green and 
Butkovic [20' , and lead to a pseudopolynomial method for finding solutions to such 
systems. This will be discussed in the last section of the paper. We also note here 
that cyclic projectors are special case of the multiplicative version of the min-max 
functions studied in [131 [T4l 141] . 

The geometry of max cones can be thought of as a special case of the multiorder 
convexity, a concept introduced by Martinez-Legaz and Singer [38j . Although this 
idea was made explicit only recently in a work by Ni^ica and Singer [39^, it is 
closely related to the set-covering conditions ioi A® x = b systems in max algebra 
[2t 1191 I45j . The multiorder principle, see Propositions 12. II and 12.31 below, leads to 
easy proofs of many statements concerning generators, extremals and bases of max 
cones, see Butkovic et al. jllj , including the tropical Caratheodory theorem, and 
Minkowski's theorem about extremals of closed cones (also Gaubert and Katz [24J). 
The multiorder principle is also important for the tropical convexity approach, 
meaning works of Develin, Sturmfels, Joswig, Yu et al. i6i, i22l[32| . since it describes 
max cones as intersections of staircases, and their extremals as elements of bases 
of monomial ideals. 

Yet another approach to the geometry of max cones, though strongly related 
to the previous one, is to represent max cones as cellular complexes, or, roughly 
speaking, as unions of ordinary convex cones. This approach was put forward by 
Develin and Sturmfels [22j , and called cellular decomposition. The atoms of this 
decomposition are well-known to specialists in convex geometry and combinatorics, 
see Joswig and Kulas [33] for more details. As it was noticed in [431, these atoms 
are column spans of uniquely defined Kleene stars, a fundamental concept in max 
algebra. 

The aim of the present paper is to bring together some geometric and algebraic 
ideas discussed above. Section 2 discusses the multiorder principle and related 
results. In particular, we show that this principle leads to a convenient new repre- 
sentation of the nonlinear projectors mentioned above. In Section 3 we recall the 
concept of Kleene stars and examine the role of their column spans called Kleene 
cones as building blocks in the Develin-Sturmfels cellular decomposition. Further 
we show that, in turn, the cellular decomposition gives rise to new max-algebraic 
objects which we call row and column Kleene stars. We relate these new con- 
cepts to the max-algebraic pseudoinverses of matrices and to tropical versions of 
the colourful Caratheodory theorem. In Section 4 we generalize the alternating 
method of Cuninghame-Green and Butkovic |20j to the case of multisided systems 
A'-^') i^ x^ — . . . = A^'^^ ® x'^. We give a proof, based on the cellular decomposition 
idea, that if the system has a positive solution, then the method converges to a 
positive solution in a finite number of steps. We also present new bounds for the 
number of iterations in the max-plus integer case, and in the general case when 
there are no solutions, in terms of the Hilbert projective distance between max 
cones. 
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2. The role of multiorder 



2.1. Generators, bases and extremals of msLx cones. Let S C M^. A 

vector u G R" is called a maa; combination of S* if 

(2.1) it = 0A„z;, A„ eM+, 

where only a finite number of A„ are nonzero. The set of all max combinations ()2.ip 
of S will be denoted by span(S'). Evidently, span(S') is a max cone. If span(S') = V, 
then we call S a set of generators for V and say that V is generated, or spanned, by 
S. In particular, the set of all max combinations of columns of a matrix A will be 
denoted by span(A) and called the column span of A. If none of the elements of a 
generating set 5 of a max cone V can be expressed as a max combination of other 
elements, then S is called a (weak) basis of V. 
A vector w G T^ is called an extremal of V, if 

w = M © w, u, w G V" => u = It or V = w. 

Extremals are analogous to extremal rays of convex cones. If v is an extremal of V 
and A > 0, then Xv is also an extremal. 

For all i = 1, . . . , n define the following preorder relation. 

u <j V <^ uuj^ < vv~^ , Uj ^ 0, Vj 7^ 0. 

The classes of proportional elements (i.e. rays) are the equivalence classes of 
these preorder relations. The importance of these relations for the geometry of max 
cones is expressed by the following principle. Denote supp(2/) := {i \ yi ^ 0}. 

Proposition 2.1. Let V — span(S'), S C R^. Then the following are equiva- 
lent. 

1- y&v. 

2. For all j G supp(y) there exists v G S such that v <j y. 

This principle appeared as a set covering condition, see Proposition l3.12l below. 
already in the works of Vorobyev [45] and Zimmermann [48, , and in the above form 
(or with a subtle difference) it appeared quite recently in the works of Joswig [32j , 
Ni^ica and Singer [39j . and Butkovic et al. [llj . see also [9] and [22] . 

As it was remarked by Ni^ica and Singer |39j . the above proposition means 
that the geometry of max cones is a special case of the multiorder convexity |38| . 
In the multiorder convexity, one has a set of order relations, and a point y is said to 
belong to the convex hull of 5, if for any order there is a point in S which precedes 
y with respect to that order. 

The following proposition is the Tropical Caratheodory Theorem, see Helbig 
[28] . Develin and Sturmfels [22], and also [HI [25]. Note that it follows from 
Proposition 12. II 

Proposition 2.2. Let 5 C K" . Then y G span(5) if and only if there exist k 
vectors v^ , . . . ,v^ € S, where k — \ supp(j/)|, such that y G span(ti"'^, . . . , v''). 

The multiorder principle also means the following description of extremals [11| . 

Proposition 2.3. Let V C R^ be a max cone generated by S and let v e V, 
V ^ 0. Then the following are equivalent. 
1. V is an extremal in V . 
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2. For some j S supp(w), v is minimal with respect to <j in V . 

3. For some j 6 supp(w), v is minimal with respect to <j in S . 

Propositions 12.11 and 12.31 lead to a number of statements about generators, 
extremals and bases of max cones [llj , we mention only the following two of them. 
An element u G M" is called scaled, if ||u|| = 1, where || • || denotes some fixed 
norm (say, the ordinary norm or the max norm) . For the following proposition see 
Butkovic et al. 11 , and also |22|, 146] for closely related statements. 



Proposition 2.4. Let E be the set of scaled extremals in a max cone V C 
and let S C R" consist of scaled elements. Then the following are equivalent. 

1. The set S generates V and none of the elements in S are redundant. 

2. S — E and S generates V . 

3. The set S is a basis for V . 



Proposition 12.41 means that if a scaled basis of a max cone exists, then it is 
unique and consists of all scaled extremals, i.e., all the elements that are minimal 
with respect to some preorder relation <i. In particular, a scaled basis of a finitely 
generated max cone V exists and is unique, and the cardinality of this basis will be 
called the max-algebraic dimension of V . 

The following result is analogous to Minkowski's theorem about extremal points 
of convex sets, and was obtained independently by Gaubert and Katz |24j and 
Butkovic et al. [11] . 

Proposition 2.5. Let V C M" be a closed max cone. Then V is generated by 
its set of extremals, and any vector in V is a max combination of no more than n 
extremals. 

Note that any finitely generated max cone is closed f| lll [32]). One may also 
think of colourful extensions of Propositions 12.21 and 12.51 in the sense of Barany |4] , 
and progress in this direction is due to Gaubert and Meunier [25] , see also Theorem 
IX^ below. 

2.2. Projectors and separation. Given a closed max cone V C M" , we can 
define a nonlinear projector Py by 

(2.2) Pv{y) ^m8ix{v eV \v <y}. 

This operator is homogeneous: Pv{Xy) — XPv{y), isotone: y^ < y^ ^ Pv{y^) < 
Pv{y^), nonincreasing: Pv{y) £ 2/, and continuous, see |16j for the proof. For 
any vector y there are coordinates which do not change under the action of the 
projector: Pv{y)i — yi- These coordinates will be called sleepers. Projectors lead 
to separation theorems of the following kind, see [161 I22L I26L 132] and introduction 
for some historical remarks. 

Theorem 2.6. Let V C M" be a closed max cone and let y £ E" be not in 
V . Then there exist a positive vector y and a max cone V ^ V containing positive 
vectors such that the set 

n n 

(2.3) H = {v\ ^y^\ > ^{Pvm7'v^} 

2=1 i=l 

contains V but not y. If y is positive and V contains positive vectors, then one can 
take y = y and V = V. 
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The set Ti defined in (j2.3|) is an instance of the max analogue of a halfspace, 
which is generally a set of the form {v \ 0"^^ u\vi > 0"^^ ufwi}. 

By comparing this to (|2.3p we see that a separating halfspace has both u^ and 
V? positive and u^ < u^, so that the inequality in ()2.3p can be replaced by equality: 

n n 

(2.4) H^{v\ ^y-'v, = ^{Py{y)Vv^}- 



The relation of Theorem 12.61 to the multiorder principle was made explicit by 
Joswig [32| . Denote, for any positive y, ^i{y) = {u G R" | u <i y}. Observe that 
[Jl^i ^i{y) — R" , and that the separating halfspace defined by (|2.3p or equivalently 
(12.41) can also be written as 



(2.5) n^ U A,iPy{y)), 

where s\{Py,y) is the set of sleepers, i.e., the indices k such that {Py{y))k = IJk- 
Thus, in terms of the multiorder, the separation theorem says that, given a point y 
and a closed max cone V, there is a point Pyiy) such that the union of some sectors 
Ai{Py{y)) contains the whole V while the complement of this union contains y. 

If a max cone is generated by the columns of a matrix A C M"^™, then, 
denoting Pa := -Pspan(A), we deduce from (|2.2p that 

(2.6) PA{y)^AcS>(A®'y), 

where A is the Cuninghame- Green inverse oi A defined by Uij — aj^ , and (g)' denotes 
the min-times matrix product. When calculating (|2.6p . we put by convention that 
0~^ = oo and 0(8'+cx) = 0. In this form (|2.6p . the nonlinear projectors were studied 
by Cuninghame- Green [19) . We also note that formula (|2.6p represents a projector 
as a min-max function in the sense of |13|, I14|, 141) . with addition being replaced 
by multiplication. 

When V is an arbitrary closed max cone, Py can be expanded in infinite sum 
of 'elementary' projectors using the following 'scalar product', or an instance of 
residuation [151 116) : 

y/v := min yiV~ = max{A | Xv < y}. 
Namely, 

(2.7) p^(y)^0y/„„. 

Formula (|2.6p is a special case of (|2.7p . when V is finitely generated. Using the 
multiorder, we can obtain the following refinement of p.7p . Denote by A the 
componentwise minimum of vectors in R" . 

Theorem 2.7. Suppose that V C M^ is a closed max cone. Then for any 
y G IR" , the components {Py(y))i, for i e supp(j/), are equal to 

(2.8) {Pv{y))^=@y/vv,, 

vGEi 

where Ei is the set of scaled points ofV, minimal with respect to <i. The projector 
Py is linear with respect to the componentwise minimum A if and only if every set 
Ei is a singleton. 
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Proof. Writing (12.71) componentwise, we have that 

{Pv{y))i = max {vi min yk.vZ^) = max min ykivkV^^y^ . 

By ProDOsition l2.5[ any closed max cone has a scaled basis E. Denote by Ei the set 
of scaled vectors minimal with respect to <i, then for all w G y and any i G supp(w) 
there is v^ G Ei such that u* <i v and hence (wfe(w*)~^)^^ > {vkv~^)~^ for all k. 
This proves p.Sp . and (|2.8p implies that if all the sets Ei consist of one element, 
then the projector is expressed by a min-times matrix. Now suppose that there is 
an i such that Ei has at least two elements, say, u and v. Then P\/(u) = u and 
Py(i') = V. If the projector is linear with respect to the componentwise minimum 
A, then Pv{uu~^ A ot~^) = uu~^ A vv~^ , hence w = uu~^ A wf ~^ S T^. As Wi = 1, 
we have that w <i v and w <i u. As u and v are both minimal with respect to <i, 
w is not equal to either of them, which leads to a contradiction with the minimality 
of u and v. The proof is complete. D 

3. The role of Kleene stars 

3.1. Kleene stars and Kleene cones. We start this section with some neces- 
sary definitions. Let A = (fly) e K'j:''". The weighted digraph D^ = {N{A),E{A)), 
whose nodes are N{A) — [n] and whose edges E{A) = N{A) x N{A) have weights 
''J^ihj) = o-iji is called the digraph associated with A. Suppose that tt = (?i, ..., ip) 
is a path in Da, then the weight of tt is defined to be w(7r, A) = ai-^i^ai^i^ . . . ai _-^i 
if p > 1, and if p = 1. A path which begins at i and ends at j will be called an 
* ~* J path. If the starting node of a path coincides with the end node then the 
path is called a cycle. 

A path TT is called positive if w(7r, A) > 0. If for all i.j e [n] there exists a 
positive i —^ j path, then A is called irreducible. 

The maximum cycle geometric mean of A, further denoted by A(v4), is defined 
by the formula 

X{A) = max //(cr, A), 

where the maximisation is taken over all cycles in the digraph and 

fiia,A)^wia,Ay/'' 

denotes the geometric mean of the cycle a = (ii, ..., i^, ii). 

The following fact was proved by Carre [12] , see also [2l I19j . 

Proposition 3.1. Let A e M"'*". The series 

(3.1) A* = I®A®A^®... 

converges to a finite limit and is equal to I(BA(B. ■ .©A"^^ if and only if X{A) < 1. 
In this case also A(^*) < 1. 

The matrix series A* defined by p.ip is called the Kleene staroi A, which comes 
from the theory of automata, see Conway |17| . Kleene stars enjoy the property 
(A*)'^ = A*, i.e., they are multiplicatively idempotent. Their diagonal entries are 
all equal to 1, i.e., the Kleene stars are increasing. Actually these two properties 
are also sufficient for a matrix to be a Kleene star, and further by a Kleene star we 
will also mean any matrix with these two properties. We also note that (A*)'^ = A* 
implies that (A*)* ^A*. 
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A max cone will be called a Kleene cone if it can be represented as max-algebraic 
column span of a Kleene star. 

In terms of the multiorder, we can say that a matrix ^ is a Kleene star if and 
only if an = 1 for all i G [n] and A.i <i A.^ for all i, k such that a^fc ^ 0. That is, A 
is a Kleene star if and only if an = 1 and A.i is the unique minimum of span(^*) 
with respect to <i for all z e [n], so that all the sets Ei defined in Theorem 1 2. 71 are 
singletons. The last sentence of Theorem 12.71 can be now formulated as follows. 

Proposition 3.2. Py is a min-times linear operator if and only if V is a 
Kleene cone. If V — span(A), where A is a Kleene star, then Pv{y) — A®' y for 
all y. 

Kleene stars play crucial role in the description of max-algebraic eigenvectors 
and subeigenvectors of nonnegative matrices. If for some x and A we have that 
A® X — Xx, then A is a max-algebraic eigenvalue of A, and a; is a max-algebraic 
eigenvector associaXed with this eigenvalue. Analogously, x is called a max-algebraic 
subeigenvector associated with A, ii A® x < Xx. 

The well-known Perron- Frobenius theorem has a max-algebraic analogue [2l [3l 
[191145] . 

Theorem 3.3. Let A e Rl'''^. 

1. A has a max-algebraic eigenvalue, and the number of such eigenvalues is 
less than or equal to n. 

2. X{A) is the largest eigenvalue of A. 

3. If A is irreducible, then X{A) is the unique max-algebraic eigenvalue of A 
and all eigenvectors associated with X{A) are positive. 

The set of eigenvectors associated with a fixed eigenvalue A is a max cone, and 
analogously the set of subeigenvectors associated with a fixed A is a max cone, 
so they will be called the eigencone and the subeigencone associated with A. For 
a nonnegative square matrix A £ M"^" the eigencone associated with 1 will be 
denoted by V{A), and the subeigencone associated with 1 will be denoted by V*{A). 
A matrix A S R!^^" is called definite, if X{A) — 1. We do not lose much generality 
when considering definite matrices, as for any matrix A with X{A) ^ 0, the matrix 
A/X{A) is definite and has the same eigenvectors and subeigenvectors as A. 

Any subeigencone is a Kleene cone, and the other way around. 

Proposition 3.4. Let A e M"''" be definite, thenV*{A) = V{A*) = span(A*). 

Proof. First note that by Proposition 13. 1[ if X{A) = 1 then A* exists and 
A(A*)^1. 

We show that V*{A) = V{A*). Suppose that A* (g) x = x, then A (® x < x, 
because A < A* . li A® x < x, then {I ® A) ®) x = x and also A* (g) x = x, since 
A"^ (E) X < x for any m (due to the isotonicity of matrix multiplication) . 

We show that V{A*) — span(^*). It is immediate that V{A*) C span(yl*), as 
V{A) C span(A) for any matrix A. If A* converges, then A igi A* ^ A® A^®..., 
so A® A* < A* meaning that each column of A* is a subeigenvector of A. Hence 
s-paii{A*) CV*{A). D 

The positivity of subeigenvectors is addressed in the following observation. 

Proposition 3.5. Let A C E"^" be such that au = 1 for all i E [n]. Then 
V* (A) contains a positive vector if and only if A is definite. 
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Proof. The "if" part: If A is definite, then by Proposition [Q V*(A) = 
span(y4.*) and we can take, for a positive subeigenvector of A, any max combination 
of all the columns of A* with positive coefficients. 

The "only if" part: Suppose that there exists a positive x such that A(^x < x, 
and take a cyclic permutation r = (ii, . . . , ifc) of a subset of [n\. Then we have that 
cLiiii^iXii^-^ < Xii for Z e [A:], assuming iu+i '■= ii- Multiplying all these inequalities 
and cancelling the coordinates of x we have that 'w{t,A) < 1. Hence X{A) < 1. As 
all diagonal entries are equal to 1, we have that X{A) — 1. D 



Proposition 13.41 implies that if A is a Kleene star, then 

span(v4) = V{A) = V{A*) = V*{A) = {x \ a^^Xj < x„ i,j G [n]}, 

and it is not hard to see the following. 

Proposition 3.6. Let K be a max cone in M" . Then it is a Kleene cone if 
and only if for some matrix B it is the solution set of the system of inequalities 
bijXj < Xi, i,j(z [n], satisfied by at least one positive x. 

Proof. The "if" part: If the system is satisfied by a positive x, then bu < 1 
for all i g [n]. Take B := I Q) B, then B has all diagonal entries equal to 1, 
K = V*{B) and there is a positive x G V*{B). By Proposition 13. 5( B is definite, 
and by Proposition l3.41 K = span((i3)*). 

The "only if" part: If iiT is a Kleene cone span(j4*), then by Proposition 13.41 
and Proposition 13 . 51 we can take B := A* . D 



The above observations imply that Kleene cones are convex cones, and that 
they have many close relatives in the realm of combinatorial geometry, see Joswig 
and Kulas [33^ . 

One may think of various systems of inequalities describing the same Kleene 
cone. However, the Kleene star which defines this cone is unique 



Proposition 3.7. Suppose that A and B are two Kleene stars. Then A = B 
if and only i/span(^) = span(i?). 

We now describe the bases of V{A) and V*{A), for a definite matrix A e M"^". 
The cycles with the cycle geometric mean equal to 1 are called critical, and the 
nodes and the edges of Da that belong to critical cycles are called critical. The 
set of critical nodes is denoted by Nc{A), the set of critical edges is denoted by 
Ec{A), and the critical digraph of A, further denoted by C{A) = {Nc{A),Ec{A)), 
is the digraph that consists of all critical nodes and critical edges of Da- All cycles 
of C{A) are critical [2]. For two vectors x and y, we write x '^ y if x = Xy for 
A > 0. The following theorem follows from well-known results on the max-algebraic 
spectral theory [2], [T9l. [23] . 

Theorem 3.8. Let A e R^^^ be definite, and let M{A) denote a set of indices 
such that for each strongly connected component of C{A) there is a unique index in 
M(A) which belongs to that component. 

1. The following statements are equivalent: A* ^ A*j, A*_ ^ A*,, i and j 
belong to the same strongly connected component of C{A). 

2. Any column of A* is a max extremal o/span(A*). 
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3. The subeigencone of A, which is the eigencone of A* , is 
V*{A) = ViA*) = J a,A* © ajA*^, a,,a, G : 

yieM(A) 3iC(A) 

and none of the columns of A* in this description are redundant. 

4. The eigencone of A is 

[ieAf(A) 
and none of the columns of A* in this description are redundant. 

Proposition 12.41 and Theorem 13.81 imply that extremals of V* [A) are precisely 
the columns of A* , so the columns of A* , after eliminating the proportional ones, 
constitute the basis of V*{A) = span(A*), and the columns whose indices belong 
to C{A) constitute the basis of V{A). Denote by nc{A) the number of strongly 
connected components in C{A), and denote by Nc{A) the set of nodes that are not 
critical. Theorem 13.81 yields the following corollary. 

Proposition 3.9. For any definite matrix A e R"^", the max-algebraic di- 
mension of the subeigencone of A is equal to nc{A) + \Nc{A)\. The max-algebraic 
dimension of the eigencone is equal to nc{A). 

Kleene cones are both convex cones and max cones. They are inhabitants of 
two worlds, that of max algebra and tropical convexity, and that of nonnegative 
linear algebra and ordinary convexity. One might think of an interplay between 
these worlds. For a definite matrix A, define the linear space 

(3.2) UC{A)) = {xe K" I a,,x, = x,, (z, j) G E,{A)}- 
A proof of the following theorem can be found in j44j . 

Theorem 3.10. Let A e R"^" be a definite matrix. Then L{C{A)) is the 
linear hull of the convex cone V*{A). The linear dimension of V*{A), i.e., the 
dimension of L{C{A)), is equal to the max-algebraic dimension ofV*{A), i.e., to 

na{A) + \Nc{A)\. 

The intersection of Kleene cones is again a Kleene cone. More precisely, we have 
the following proposition, see Butkovic [7J for the case k — 2. The proof is based on 
the formula (A* ©S*)* = (A* (g)B*)*, which follows from (A©S)* = A* (g) {B (g) A*)* 
[17j . and on the observations above. 

Proposition 3.11. Let A'^^\. . . , A'^''^ e R"""" be Kleene stars. The following 
are equivalent. 

1. pl^"^-^ span(yl('') contains a positive vector. 

2. A(etiA(')) = i. 

3. A({^j^j^ ^('^(*))) — 1 for some permutation tt of {1, . . . , fc}. 

4. A({g)^^i ^(^(«))) ^ 1 for all permutations tt of {I, . . . ,k} . 
If any of these equivalent conditions are true, then 

k k k 

(3.3) fl span(A«) = span((0 ^W)*) = span(((g) ^(-W))*) 

i— 1 2—1 i— 1 
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for all permutations ir. 

Proof. Complete Kmax, x with +00 and assume a x +00 = +c» for any positive 
a and x +cxd = 0. Matrix algebra over this completed semiring is a regular algebra 
in the sense of |17| . This means in particular that A* is always defined, (A*)* = A*, 
{A © B)* = A* (g>{B(g) A*)* and {A (g) B)* = I ® {A (g) {B (g) A)*). If A and B are 
two Kleene stars, then 

(AgiB)* =I®{A(g,{B(g, A)*) ^A(g){B(g,A)*^ 

^^■^^ ={A®By = {B®Ay = {Brg,A)*. 

It can be shown by induction that (A^^)®. . .©AC'))* = {A'^^^'^^g). . .g)A'^''^^'>)* for any 
permutation tt of {1, . . . , /c}. Using Proposition l3. ll we obtain that A(0j^^ A^*' ) < 1 
is true if and only if A({^j^j^ A'^^^^) < 1 is true for some tt, and hence if and only 
if the same is true for all tt. The inequalities here can be replaced by equalities, 
since all diagonal entries, and hence all eigenvalues, of any product or entrywise 
maximum of Kleene stars, are greater than or equal to 1. This yields equivalence 
of 2., 3., and 4. 

We now prove the equivalence between 1. and 2., and (|3.3p . We have that 



(3.5) y*(0 A«) - fl y*(^«) = fl span(AW), 

i—l i—1 i—1 

where the first equality is immediate, and the second equality follows from Proposi- 
tion [Jill Note that all diagonal entries of ©j^^ ^^*-' are 1, and by Proposition l3.51 
^*(®i=i A^^^) contains a positive vector if and only if A(0,-^-^ A^*') = 1. This, to- 
gether with (13. 5p . implies the equivalence between assertions 1. and 2. By Propo- 
sition [33 l/*(0ti^^'^) = span((0ti^^'^)*) since X{®ti^^'^) = 1- which 
yields dSSl). □ 

3.2. Cellular decomposition. We have described some properties of Kleene 
cones. Though such cones are very special, they can be viewed as building blocks, 
or atoms, of any finitely generated max cone. This can be seen as the main idea 
of the cellular decomposition, an ingenuous concept of Develin and Sturmfels [22j . 
which we adjust below to the setting of max cones. 

Let A C M"^™ be a nonnegative matrix with m nonzero columns and n nonzero 
rows. The column type of y with respect to A is defined to be the ?7i-tuple of subsets 
Ti , . . . , Tm of [n] , where every Tj , for j £ [m] is defined by 

Tj ^{ie [n] I fty-yri > akjy'^^, k e [n]} ^ {i e [n] \ y >i A.^}. 

The row type of y with respect to A is an n-tuple of subsets Si, . . . , Sn of [m] , where 
every Si, for i € [n], is defined by 

Si = {j e H I Oijy^r'^ > akjy^^, k e [n]} = {j e [m] | y >, A.j} = 

The theory oiAg)x^y systems [H [H [l9l [22l HH |48] is based on the following set 
covering conditions for y to be in span(^), see the proposition below. The multi- 
order principle fProposition l2.ip can be seen as a reformulation of these conditions, 
therefore we leave the proposition below without proof. 
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Proposition 3.12. Let A e M"^™ have all rows and columns nonzero and let 
y G M" be a positive vector with the column type T ~ (Ti, . . . , Tm) and the row type 
S — (S*!, . . . , Sn)- The following are equivalent. 

1. y e span(yl); 

2. [jtlT^ = [n]; 

3. none of Si,i £ [n] are empty. 

See also Akian et al. [I] for an infinite-dimensional generalisation in the context 
of Galois connections. 

Following Develin and Sturmfels [22j . we can see this from a geometric view- 
point. For any row type S, we define its region with respect to A by 

Xs = {y positive | yky~^ > akja^/, yk,i, Vj £ Si}. 

Proposition l3. 121 means that the part of span(A) consisting of all positive vectors is 
the union of the regions Xs such that S do not contain empty sets ( [22j . Theorem 
15). If Xs is not empty, then the closure of Xs is 

(3.6) cl{Xs) ^{y&Rl\ akjoT^^y, < y^, Vfc, i, Vj G S,}. 

It follows from the results of '22' that the relative interiors of regions build up a 
cellular decomposition of the positive part of W^. We will need a weaker statement, 
but without positivity. 

Proposition 3.13. Suppose that A e R"^" has all rows and columns nonzero. 
Then the max cone span(A) is the union of c\{Xs) such that Xs are not empty and 
S do not contain empty sets. 

Proof. As A has all rows nonzero, the max cone span(yl) contains positive 
vectors. By Proposition l3 . 1 21 if y is positive, then y G span(A) if and only if the row 
type of y does not contain empty sets. Hence the positive part of span(A) is the 
union of nonempty Xs such that S do not contain empty sets. Further, span{A) is 
the closure of its positive part. Indeed, span(>l) contains positive vectors and for 
any u G span(^) and a positive v G span(A) we can take w = u (B sv G span(A) , 
so that llw — m|| < e\\v\\ (the max norm) and w is positive. Hence span(A) is the 
union of closed regions cl{Xs) such that Xs are not empty and S do not contain 
empty sets. D 

From the max-algebraic point of view, an important role in the cellular decom- 
position is played by strongly definite matrices, which are definite matrices with all 
diagonal entries equal to 1. Note that any Kleene star is a strongly definite matrix. 

Observe that cl{Xs) is the subeigencone of the nxn matrix A^ = (afj) defined 

by 

^s _ /©fees^a^fc^fe^' if ^j 7^ 0, 
Sij, if Sj = 9, 

where Sij are Kronecker symbols {Sij — ii i ^ j and Sij ~ I if i — j). It is 
immediate that all diagonal entries of A^ are equal to 1. We have the following 
proposition which can be used to compute the generators of any closed region, a 
preliminary version of this proposition appeared in [43] . 

Proposition 3.14. The closed region cl{Xs) contains positive vectors if and 
only if A is a strongly definite matrix, and in this case cl(X5) = V*{A ) = 
span((A^)*). 



(3.7) af^ 
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Proof. From ([3^ and (lajjl one infers that d{Xs) == V*{A^). After that, 
the claim foUows from Proposition 13.41 and Proposition 13.51 D 



Propositions 13.131 and 13.1^ have the following consequences. 

Proposition 3.15. For any matrix A G M"^™ with no zero rows there exist 
Kleene stars A^^\ . . . ,yl(') e IR'|''" sueh that span(A) = ULi span(A(^)). 

Proposition 3.16. For any matrix A S M" ™ with no zero rows there exist 
Kleene stars A'-'^\ . . . , A'-^^ e M"^" such that for any y e R^ we have that PaV = 
^(fc) (g)' y for some k. 

To express the dimension of a region, Develin and Sturmfels I22 j introduce the 
undirected graph Gs'- The set of nodes of this graph is [n], it contains all loops 
(j,i), and for i ^ j an edge (i,j) belongs to Gs if and only if there exists /c G SiHSj. 
The following observation relates this notion to max algebra. 

Proposition 3.17. Let A G K."^™ be a matrix with no zero rows and columns, 
let y G M" be a positive vector and S be the row type of y with respect to A. Then 
Gs^CiAS). 

Proof. Note that as all entries of A^ are equal to 1, the graph G{A^') contains 
all loops. 

Let i ^ j and (i,j) G Gs^ then there exists fc G 5; n Sj. It follows that 
atkaj,} = yiy~^ > auaji^ for all I G Sj, and therefore a^- = aikaj,}. Analogously, 
a^j = ajka~^, and therefore af^a^^ = 1 so that (i, j) G G{A^). 

Let (i, j) G C{A^), then observe that a^y^ < yi is impossible, because the 
multiplication with other inequalities over the critical cycle would lead to 1 < 1. 
So a^-2/j — yi, and hence there exists k G Sj such that Oikajj^ yj = yi. But then 
also k £ Si and {i,j) G Gs- □ 

The equality Gs — G{A^) means that G{A^) is symmetrical and {i,j) G Gs if 
and only if {i,j) or equivalently (j, j) belong to G{A^). Theorem 13.101 and Proposi- 
tion [OT] yield the following result, see also Develin and Sturmfels [22j . Proposition 
17. 

Theorem 3.18. Let A G R"^™ &e a matrix with no zero rows and columns, 
let y be a positive vector and S be the row type of y with respect to A, then both 
max-algebraic and linear dimensions ofc\{Xs) are equal to the number of connected 
components in Gs- 

3.3. Row and column Kleene stars. For a matrix A — (aij) G M"^" and 
any permutation a £ Sn (where Sn denotes the group of all permutations of [n]) 
define the weight of a to be w{a) :— Y[i=i ^iaU)- The max-algebraic permanent of 
A is defined as 

(3.8) per(A) = u;(a), 

and a permutation, at which the maximum in (|3.8p is attained, is called a maximal 
permutation. For any permutation a, define the diagonal matrix D°' = (d^) by 

if i = cr{i); 
otherwise. 
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Observe that A{D'^)^^ is an instance oi A^ , for the type S — {{cr(l)}, . . . , {a{n)}}. 
The subeigencone V*{A{D'')^^) is precisely the closed region c\{X^). It contains 
positive vectors if and only if A{D'^)^^ is strongly definite, and this is true if and 
only if the permutation a is maximal [9 . This is also equivalent to {D'^)~^A being 
strongly definite. Further A{D°')~^ will be denoted by A'^"' and {D°')^^A will be 
denoted by A^ 



(3.9) <J = a,^(j)a^.^(^.), </ - a<,-i(,),aff-i(i)j- 

The Kleene stars of A"'' and A"" will be denoted by A"''* and A'""* and called 
column Kleene stars and row Kleene stars, respectively. 

The results of Yoeli |47j . see also Cuninghame- Green J19j . Theorem 27-11, and 
Izhakian [291, I31| suggest that row and column Kleene stars are related to the max- 
algebraic pseudoinverses of matrices. The pseudoinverse of A is defined, see |19| 
and [47], as A^ = (per(^))~M'*'^J. Here A^'^^ is the pseudoadjugate of A defined 
by a^ ■" = per(v4ji), where Aji is the complementary minor to Uij. The following 
proposition collects some facts about strongly definite matrices, which are due to 
Yoeli and Cuninghame-Green. 

Proposition 3.19. Let A e K"^" be strongly definite. 

1. I <A<A^ < ... < A"-i = A" = . . .. 

2. A* = A"-i. 

3. A* = A^dj =^v. 



Izhakian |29l I31j studies the products A (g) A^ and A^ (g) yl over extended 
tropical semiring, with the main emphasis on the questions of regularity and rank. 
In this context, he proves [31] that the products A (g) A'^ and A"^ (g) A are Kleene 
stars. Below we give an elementary proof that over max algebra, these products 
are equal to column and row Kleene stars, respectively. 

Theorem 3.20. Let A G M"^" have nonzero permanent. For any permutation 
a with maximal weight we have that A""* = D'^A^ == A'SiA'^ and A"'* = A^D'^ = 
A^ ®A. 

Proof. Using (|3.9p and the definition of A^'^\ we write: 

%■'' ^ VF li"feir(/c)= (p 11 «CT-i7r(fc),7r(fe)afe'^cr-i7r(fe) = 
7r:7r(j)=i k=^j TT:7r{j)=i k^j 

fc#i 7r:7r(j)=i k^j 

= per(A) . a;_\(,), • [] «r.(.) = Per(^) • «;-\(.).(«^^)fA(,),. 

7r:7r(j)=cr-i(i) k^j 

By Proposition [XTOl (A'='^)adj ^ j^ca* ^ ^^ ^^ j^^^^g obtained that A^'^J = 
per(yl)(D'^)-iA^'"*, and hence A^ = (£i'")-i(A='")* and D^A"^ = {A"")*. We 
now infer that 
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Thus A ® A^ = A'^'^* . On the other hand, one can similarly obtain that 



^'^i — T^,=.^( A\r,i'(y* ^-1 



a-- = per(A)a[_^(^)a7;(^.) and that A^ ®A^ A^ D'' = A'-''*. D 

Clearly this theorem yields the following corollary the first part of which was 
obtained in |43) . This corollary means that for any matrix with nonzero permanent, 
both row Kleene star and column Kleene star are uniquely defined. 

Proposition 3.21. Let A G R"^" have nonzero permanent. Then for all 
permutations a with maximal weight, the corresponding column Kleene stars A'^"* 
are equal to each other, and the row Kleene stars A^'^* are also equal to each other. 

The idea of the proof in |43j was to notice that the (sub)eigencones of A'^" 
are the same for all maximal permutations a, and to use Proposition 13 . 71 that any 
Kleene star is uniquely defined by its column span. 

For a square matrix A, the span of its column Kleene star is the only region 
of span(j4) which may have full linear dimension, and the linear dimension of that 
region determines the tropical rank of A, introduced by Develin et al. [21) . and also 
investigated by Izhakian [301 . When the tropical rank is full, the interior of span of 
the column Kleene star is the simple image set of A studied by Butkovic [8| : It is 
the set of vectors y G M" such that Ax — y has a unique solution. In what follows, 
the span of column Kleene star of A will be called the essential span of A. 

The following theorem, which is a slight generalization of Theorem 8 by Gaubert 
and Meunier |25| , illustrates the role of essential span in the geometry of max cones. 
It can be thought of as a colourful generalization of Minkowski's theorem for max 
cones in the sense of Barany [4 . 

Theorem 3.22. Let U CWl be a closed max cone and let V'^,...,V" C R!^ 
be closed max cones such that the intersection of V^ with U is nontrivial for all 
i e [n]. Then there exist vectors v^ , . . . , u" such that v^ is an extremal of V^ , for 
i G [n], and span(u"'^, . . . ,v") has nontrivial intersection with U. 

Proof. Take any nonzero points y^ G F^ n C/, ...,?/" G ¥'"■ D U and consider 
the matrix A G R"^" with columns A.i = y*, for i — 1, . . . ,n. Assume first that 
A has permutations with nonzero weight. The essential span of A is the closed 
region cl(X5), where S = {{cr(l)}, . . . , {<T(n}}, for any maximal permutation a. 
Take any u G cl{Xs), then u G U and u >i ^.^.(i) for all i. The column ^.^(i) is 
equal to y'^^'' and it belongs to V^^"^' . Applying Minkowski theorem (Proposition 
12. 5|) and the multiorder principle (Proposition [23| , we obtain an extremal i''^'*' of 
y"'^') such that w'^'^*^ <i y'^'^*^ <i u. Applying Proposition 12.11 again, we see that 
u G span(t;'^(^\ . . . , w'^(")). As u G f7, the claim follows. 

In the case when A does not have nonzero permutations, an inductive argument 
using Hall's marriage theorem, see [25j . shows that there exist subsets of indices M, 
A^i and N2 such that the submatrix Ajatj^m] is zero, while the submatrix At^jsj^M] 
is square and has a permutation with nonzero weight. Then the above argument 
goes with the essential span of that submatrix. D 

4. Cyclic projectors and the alternating method 

4.1. Cyclic projectors and separation of several maix cones. Let 

V^ , . . . , V^ be closed max cones in R!Ji and denote by Pi the projector onto V^ . The 
composition Pk ■ ■ ■ Pi will be called the cyclic projector associated with V^, . . . ,V''. 
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This operator inherits many properties of the sole projector: it is a homogeneous, 
continuous, isotone and nonincreasing operator. In general, it is not linear with 
respect to max and min operations. Such operators can be treated by nonlinear 
Perron-Frobenius theory. In particular, the following theorem of Nussbaum |40| 
generalizes the well-known CoUatz-Wielandt formula for the spectral radius of a 
nonnegative matrix. 

Theorem 4.1. Let F be a continuous, homogeneous and isotone operator in 
M" . Then the spectral radius of F is equal to 

(4.1) r{F) = inf{A | 3y positive: Fy < Ay}. 

Such operators have no more than one eigenvalue over any set of vectors with 
the same support, and therefore the total number of their eigenvalues is finite. For- 
mula (|4.ip implies that the spectral radius is monotone. Define the cyclic projective 
distance oi y^, . . . ,y'' G M" by 

(4.2) pH(2/\...,y'=)-log yl{y^J-'-----yt{yl)-\ 

ii.....Ak&M 

when supp(y^) = . . . = supp(y'^) — M, and by -f-oo otherwise. In the case k = 2 
this is the Hilbert projective distance between two points in K". An equivalent 
definition is 

k 

(4.3) pH{y\ ■■■,y'') = loginflj] A, I y' < W+\i e [k]}, 

i=l 

where y''^^ := y^. Note that pn is stable under multiplication of the arguments by 
nonzero scalars and under their cyclic permutation. If Ym=i Vi ~ ^ ^'^^ ^ ^ i^]' then 
it follows from (|4.3p that Ai > 1, and ph(j/^, • • ■ , y'') = if and only if j/^ = . . . = y'^. 
For general y^, . . . ,y'^ e R" \{0}, pii{y^, . . . , y''') = if and only if y^, . . . , y'"' are 
proportional to each other. 

Define the cyclic projective distance between closed max cones V^, . . . , V^ by 

(4.4) pn{V\ ...,V'')^ ^ inf pH(y\ . . . y'). 



The minimum in (j4.4p is attained since pu is lower semicontinuous, see Proposi- 
tion |4i8] below. 

The monotonicity of spectral radius is crucial for the following theorem 



Theorem 4.2. Let V^, . . . ,V'' be closed max cones in R" . Suppose that y° is 
an eigenvector of Pk ■ ■ ■ Pi associated with the spectral radius, and consider vectors 
y^ e y\ . . . , y*-' e V'^ defined by y^ := Piy°, . .. ,y^ := Pky^~^ . Then 

PH(y\...,y') = PH(V'\...,V^')--logr(Pfc...Pi). 

Cyclic projectors also enable to prove a separation theorem for closed max cones 
[26j . with the following ideas in mind. Firstly, formula (|4.ip implies the existence 
of a positive subeigenvector with A < 1. Secondly, if we take such a subeigenvector, 
then its projections onto V^, . . . ,V'^ define separating halfspaces, see Theorem l2.6l 



Theorem 4.3. LetV^, ■ ■ ■ jV' C M^ be closed max cones. If each ofV^ , . . . ,V^ 
has a positive vector, then the following are equivalent. 

1. There exists a positive vector y and A < 1; Pk ■ ■ ■ Piy < ^y . 
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2. There exist halfspaces Ti} , . . . .H^ such that V^ C H'^, . . . ,¥>" C H'' and 

ntiH* = {o}. 

4. r(Pfe---Fi)<l. 
The statements 2. and 3. are equivalent even if V^, . . . ,V'' do not have positive 
vectors. 

4.2. The alternating method and its convergence. In what follows we 
consider the case when V^ = spa.n{A^^^), . . . ,V'' = span(^('')), and A''^\ . . . , A'-''^ 
are nonnegative matrices with an equal number of nonzero rows. A natural question 
is to find a positive solution to the system of equations 

(4.5) yl(i)«)xi = ... = ylW®x^ 

and the cyclic projectors provide an efficient method for doing this. 



ALTERNATING METHOD 

Input: Nonnegative matrices A^'^'> £ R"''™\ . . . , A^*^) e r'J'X™'. ^jth an equal 
number n of nonzero rows. 

Initialization: Arbitrary positive y*^"^ :— y^^^^^. 

Iteration: Number / > 1. For all s = 1, . . . , fc compute a;'-'-''* := A(*) ®' y(')'*"^ and 
y(i)s ._ ^(s) ^ ^(0^^ Set xW := a;^')'^ and yW := y^)*^. 

Stop: If y(') = y^'^^\ then stop. The vectors x^'-''*, for s = 1, . . . ,k, give a solution 
to system (j4.5p . Else if y\ < y\ for all i g [n], then stop. There is no solution. 



Over the semiring Mijiax,+ — (K. U {— oo}, = max, (g) = +) and for fc = 2, this 
method was formulated by Cuninghame-Green and Butkovic }20j . The method 
is essentially a max-algebraic version of the cyclic projections method known in 
optimization theory [5j, since y^'-* — Pk ■ ■ ■ J^iy*-'"^-*- 

The first part of the stop condition follows from the fact that Pi, . . . ,Pk are 
nonincreasing projectors onto span(A(^)), . . . , span(yl('^) ) . Indeed, if y*-'^^' = y*-'', 
then the inequalities 

y(')>Pfe-i---Piy('-^)>...>Piy('-^)>y('-^) 

are satisfied with equalities, implying that y^'-** = Pg ■ ■ ■ Piy'^^^ are equal for all 
s e [k] and that yC) € span(yl(i)) n . . . n span(yl('=)). As yC)" = A^''^ (g) a;(')^ for 
s £ [k], we have that x*-'^^, for s € [k], give a solution to (|4.5p . 

Also note that the absence of zero rows in the matrices implies that all vectors 
in the sequence generated by the alternating method are positive and hence any 
solution, which the alternating method may find, has to be positive. 

The following proposition, similar to the results of |20j . justifies the second part 
of the stop condition. It emphasizes the role of sleepers, i.e., such indices i{s) £ [n] 
(for 8 = I,. ..,k) that y| ?" = ylj^ = ... for the whole sequence {y^'^'', I > 1}, 

and j{s) £ [rus] such that x^J'^ — x^J'' = . . . for the whole sequence {x^'^", I > 1}. 
Sleepers will be called eternal, if the corresponding coordinates are constant for all 
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Z > 1, and temporary, if the corresponding coordinates are constant up to the last 
iteration of the alternating method. 

Proposition 4.4. Let A^^\ . . . ,A^''\ x^'-^^ and y^')"* be as in the formulation 
of the alternating method. Then 

1. temporary sleepers exist for all sequences {x'"*} and {y^'''"}, s G [k]. 

2. if (14. 5p has a solution, then eternal sleepers exist for all sequences {a;*-"*} 
and {y(')"}, s G [fc]. 

3. if (j4.5[) has a positive solution, then {x^'''''} and {y^'''''}, for all s G [k], 
are bounded from below by positive vectors. 

Proof. 1. Assume that for some s G [k] and Z > 1 we have that all coordinates 
of y^'^" or x^'^* are strictly less than that of y*-^^'* or x*^^^^. Then we have that 
yWs < ^y(i)* or a;^'-''* < pLx'^^^" for some /i < 1. As all matrix multiplications are 
homogeneous and isotone, we have that y''^ < /j.j/'-^-' so that all coordinates of y^'^ 
are strictly less than that of y'^^ and the alternating method immediately stops. 

2. and 3. Take any s G [k]. If there is a vector y in the intersection of column 
spans, we can scale it so that y < j/^^^** and yi = yi for some i. In terms of the 
multiorder, y <i y'-^-''' (for this scaling it is essential that y'^^'^ and hence y^^^* are 
positive). As the projectors are all isotone and y is their fixed point, we have that 
y < y^'-*^ and yi = yl for the whole sequence. If (|4.5p has a positive solution, then 
the same scaling argument shows that the sequence {y^^^'', j/'^^^'*, . . .} is bounded 
from below by a positive vector. Now note that the same line of argument applies 
to {a;(')^} as well. D 

In what follows we will prove that the alternating method converges to a positive 
solution if a positive solution exists. We note here that a cyclic projector is a 
min-max function in the sense of [13L 1141 141] . with addition being replaced by 
multiplication, and the convergence of the alternating method follows from the 
results of [131 141) concerning the ultimate periodicity of min-max functions. Below 
we give a different proof which uses the cellular decomposition idea. 

We first investigate the convergence of the alternating method for Kleene stars, 
which then enables us, using cellular decomposition, to prove the finiteness results 
for general matrices. 

Proposition 4.5. Suppose that A^^\...,A^^^ e M"^" are Kleene stars. If 
span(yl^^^) n . . . nspan(A^'^') contains a positive vector, then the alternating method 
converges in no more than n iterations. 

Proof. The alternating method starts with an arbitrary positive initial vector 
y and repeatedly applies the composition Pk ■ • • Pi. Due to Proposition IS . 21 we have 
that 

Pfc . . . Fiy = ]4W ®' . . . ®' ]4W ®' y, 

and hence 

{Pk--- Pirv - (^4W ®' . . . ®' IW)™ ®' y. 

This means that the stabilization of the alternating method is equivalent to the 
stabilization of {A^-^'' ... A^'^^)™ (g) y for any positive y. Denote the matrix 
product A(i) ® . . . ® AC"') by C. By Proposition [3lTl we have that A(C) = 1. We 
also have that the diagonal entries of C are equal to 1 and hence it is a strongly 
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definite matrix. By Proposition 13.191 the powers of C stabilize in no more than 
n — 1 steps, and this proves the claim. D 

Now we make use of the cellular decomposition to prove that if there is a 
positive solution, then the alternating method finds a positive solution in a finite 
number of steps. First we prove the following technical proposition. 

Proposition 4.6. Suppose that A^^\ . . . , A^^^ e M"^" have all diagonal entries 
equal to 1 and suppose that any product D of no more than n of them has X{D) < 
1. Fix a mapping j : {1,...} ^-> {1, . . . ,k}. Consider the sequence of products 
C(™) = ^0-(™)) . . . ® A(^'(i», /or m > 1. Then there exists m < n^ ~ 1 such that 

Proof. For the case of just one matrix, this is Proposition 13.191 We argue 
by induction, assuming the result is true for fc — 1 matrices and proving it for k. 
Choose any mapping tt : {1, . . . , n} i— > {1, . . . , k}. Then either for some m < n^ we 
have that there are no repetitions before that m and 

n 

4=1 

where each B''^^ is a product of less than n^^^ — 1 matrices, or there is a repetition, 
and in this case we are done. Hence, for M = n^ — 1, either there are repetitions 
before that M, or the product C*^*^^ — (q ) contains all the above mentioned 
products. We claim then that 



(4.6) 4")= a 



(m) _ ^^ jT^in)) _ ^ (^(1)) 

iin-i ■ ■ ■ n i 



for all m > M. Indeed, c\™ is greater than or equal to the maximum on the r.h.s. 
due to the choice of M and since all diagonal entries of all matrices are 1. It is 
actually equal to this maximum because all products of no more than n matrices 
have A < 1, so the weight of any path of length M does not exceed the weight of 
the simple path obtained after cycle deletion, and the weights of all simple paths 
are already in (|4.6p . D 

Theorem 4.7. Suppose that A^^^ G R+^'^S . . . , ^C^') e R'jx™'" have all rows 
nonzero and are such that span(yl'"'^') n . . . n span(A'^) contains a positive vector. 
Then the alternating method stabilizes in a finite number of steps. 

Proof. It follows from Proposition 13. 151 that for each matrix A'^'^ we have a 
Kleene decomposition 

span(A(*))== |Jspan(A(*')), 
1=1 
where A^^'-^ G R"^" are Kleene stars. Then we have that 



(4.7) 



(Pfc • • • Pi)"y = (^(feK*,™)) 0' . . . ®' ^("(1,™))) ®' . . . 



(g,' (^(fcKfc.i)) ®' . . . (g)' ^(iKi.i))) ®' y 
for some index mappings l{i,j). 
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It suffices to prove the stabilization of the sequence 

(4.8) S(™) ® . . . ® b(i) ® y, 

wliere B^^) = (-(-^(feK*:,*)))^ (g, . . . ^ (A(i'(i^*»)^). Note that the number of matrices 
B^^^ is also finite. Since the spans of the matrices A^^\ . . . ,A^''^ have a point in 
intersection, by Proposition 14.41 sequence (|4.7p is bounded from below, and hence 
is bounded from above. 



Consider a finite product B of some matrices B^-^\ appearing in (|4.8p . If X{B) > 1, 
then at least one of the matrices making this product will appear only a finite num- 
ber of times. Otherwise the sequence will be unbounded, which is a contradiction. 

Hence after some finite m the matrices B^^' appearing in the sequence will be such 
that X{B) < 1 for any product B of no more than n of them. 

After that, the finite convergence of alternating method is guaranteed by Proposi- 
tion SH D 



4.3. Bounds on the number of iterations. Now we examine the case when 
the system has no solution, i.e., when the max cones span(A'^^)), . . . , span(A('^)) do 
not have nontrivial intersection. Here we will need the total projective distance 
between y^, . . . ,y , which is the sum of projective distances 

(4.9) psiy\ •••,/) = Puiy\y^) + ■■■ + PH(^/^2/'), 

if 2/^, . . . , y*^ have equal supports, and +cxd otherwise. Note that 

(4.10) p^iy\ ...,y') = pu{y\ ...,/) + Pu{y\ ■■■, y'), 

where pn is the cyclic projective distance defined by (|4.2p . By analogy with (|4.3p . 

fc 

(4.11) p^{y\ ...,y'')^ loginf{[] Kp, \ f < X,f+\ y'+' < p,y\ i G [fc]}. 



,fe+l 



where y :— y . Like pn, the total projective distance is stable under scalar 
multiplication of the arguments and their cyclic permutation. 



Denote iS„ := {x S M" | X)"^! ^i ~ ^i ^"^^ consider S^ := 5„ x . . . x 5„ 
endowed with product topology. A function : 5^ i— > R4. U {+cxd} is called lower 
semicontinuous if the sublevel sets 

(4.12) 5,';(<^,a) = {{y\ ...,/') e S^;, \ ^{y\ . . . ,y'=) < a}, 

are closed for all a £ M.^ . The author gratefully acknowledges the idea of the proof 
of the following proposition to Stephane Gaubert. 

Proposition 4.8. p-s:{y^, ■ ■ ■ ,2/*^) and pii{y^ , . . . ^y^) are lower semicontinuous 
onS^. 

Proof. Consider sequences {y*•'"•'^ rn > \\ <Z Sn converging to y\ for i G [fc]. 
We need to show that if (y(")\ . . . , ?;("''=) e5^(ps, a) (resp. if (2/(™)\ . . . , y'^''^) G 
'5^(PH,a)) for all m > 1, then (yi,...,yfc) G 5,^(ps,a) (resp. (yi,...,y'=) G 
5^(pH,a)). If (y(™^\ . . . ,2/^"^'') e 5^(ps,a) for all m, there exist A,^'"\//|'"^ G M+ 
such that y(™)* < X^^'^ y(■m)^+l ^nd y(™)^+i < pf^y'^'^)^ for all i G [fc], and that 
nti Al"Vi"^ < a. As Er=i y^^' = 1 for all m and z, and y^™)' < X^[^\i--)-+^ 
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and 2/(™)»+i < ^f')y(™)% we have that a|™' > 1 and ^|"^ > 1. Usmg these in- 

equahties and 0^=1 '^i Vi ^ «: ^'^ obtain that 1 < Xl™' < a and 1 < fj,^' < a 
for all i £ [k]. Taking convergent subsequences if necessary, we can assume that 

^i ~^ ^i ^^^ l^i ~* Mi for * G [^]- Then we have y* < A^y'+i and z/*+^ < /iij/' 
for all i e [fc], and ni=i '^iMi 1^ c-i which yields (y^, . . . , y^) G S!^{ps, a). The proof 
for the case of ps is complete, the case of pn is treated analogously. D 

By analogy with (|4.4p . the total projective distance between closed max cones 
V^,. . . ,V^ is defined by 

p^{V\ ...,V^)= pn{V\ V^) + ... + pn(V\ V') = 

(4.13) . / 1 fe^ 

= mm pY,{y ,...,y ). 

y^ev^ ,...,y''eV'' 

Observe that ps{y^, ■ ■ ■ , y^) = if and only if y^, . . . , j/*^ are multiples of each other. 
This is generalised in the following proposition. 

Proposition 4.9. Let V^,...,V^ C Wi^ he dosed max cones. Then 
Py:{V^ , . . . ,V^) — (equivalently, pB_{V^, . . . ,V'') = 0) if and only if the inter- 
section of V^ , . . . , V^ is nontrivial. 

Proof. We show the "only if" part. The intersections of V^ and 5„ are closed 
sets. Let the sequences {j/'-™-'% m > 1}, for i G [k] and j/(™)* e V^CiSn, be such that 
lim™_oo /9e(2/(™)\ . . . , y^™)'^) - (or lim™_oo Pii{y'^"'^\ ■ ■_, y'^^'^) = 0). As <S„ is 
compact, we can assume that y(™)* -^ y* for i G [k], where y' G V^OSn as F'n5„ is 
closed. Proposition 14.81 implies that ps{y^, . . . .y^) = (resp. ph(?/^, ■ • • ,2/*^) = 0). 
Hence y', for i G [A:] , are proportional vectors contained in 1^^ n ... n V*' . The proof 
of the "only if" part is complete. The "if" part is obvious. D 

Let vector y and matrix A have finite entries. Denote 

(4.14) I |y| I =log0 y,y7i, ||A|| = log a.^aT^. 

A vector y — Ai=i ^i-^ij where A^ > for all i G [n], and A denotes the componen- 
twise minimum, will be called a min combination of the columns of A. 

Proposition 4.10. Let A G M"^" and y G M" have all entries positive. If y 
is a max combination or a min combination of the columns of A, then \\y\\ < 



Proof. Let y = 0^- XjA.j, or let y = A^ AjA.j with aU A^ ^ 0. Then 
exp(||y| I) = 0y.y-^ = 0(0 Afea.fe).( /\ X^'aJ,') ^ 

= Afefl,,. • ( f\ Xi\ji^) < a^fcaT^ < cxp(PII), or 

l,j,k l:\iT^O lJ,fc:Afe^O 

exp(||y||) ^^y.y-' - ^{/\^ka.u) ■ i^K'^j') - 

i,j ij k I 

= V'«7/' • (A^fe"^*^-) < 0«''«7/' < cxp(Pii), 

i.j,l k ijd 

respectively. The claim follows by the monotonicity of the logarithm. D 
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Proposition 4.11. Let u e M" be a positive vector, let V '^M.\ be a closed 
max cone and let v = Pv{u). Then X]j=i(l'-'S'"i ~ logfj) > Ph(u,w). 

Proof. As v < u and u^ = Vk for some k, we have that pYi{u,v) = 
max"^]^(logUi — \ogVi). As any sum of nonnegative numbers is greater than or 
equal to any of its terms, the claim follows. 

D 

Proposition 4.12. Suppose that A e R"^™, and suppose that x^,x^ E M™ 
and y^,y^ G M" are positive and such that y^ > y^ with strict inequalities in at 
most n' coordinates, x^ > x^ and A® x^ = y^ , ^ (8> a;^ = y'^ ■ Then 

1. there exists k such that x^ix^)^^ > maxj yl{y'^)~'^ ; 

2. the inequality YllLiC^ogxl - logx|) > ^ Y^-^iilogyj - logyf) holds. 

Proof. Let t be such that max^ y].{y1)^^ = yKy^)'^ and define k such that 
moKs^atsxl) = atkx\ = y}- The inequalities atk 7^ and atkx\ < y^ imply part 1. 
To obtain part 2. we recall that any sum of nonnegative numbers is greater than 
or equal to any of its terms, and that the maximum is always greater than or equal 
to the arithmetic mean. D 

Now we obtain a bound for the number of iterations of the alternating method. 
For brevity, we denote p^{A'-'^\ . . . ,A'-'''>) := /os(span(^(^)), . . . ,span(A(''))). 

Theorem 4.13. Suppose that A(i) e E'j^^S . . .,A<-''^ e ]R'|^"^ that A^^) has 
all entries positive, and that span(yl''^)) fl . . . fl span(^('^') = {0}. Then after no 
more than 

(4.15) 2(n - l)min(pW||, (m^ - l)pW^||)/ps(^W, . . . ^y^W) 

iterations the alternating method will terminate. 

Proof. Let the sequences {y^'^", I > 1} and {a;^'^'*, I > 1}, for s G [k], be as in 
the formulation of the alternating method. Using Proposition 14. 11[ we obtain the 
following lower bound for the total sum of logarithmic coordinate losses of y*-'^ at 
each iteration: 



n k—1 



(4.16) tt 



(logyp^' -logyf)) = J2J2(^gyl^'+'~\ogyt) > 



s=0 ? 



>ps(y«\...,y«'=)>p^(AW,...,AW). 
Using Proposition 14. 12[ we also obtain that 

J2(logxt'' -logxf^) > ^E(logyf+^) -logyf)) > 
(4.17) »=i " '=1 

>^ps(aW,...,a(^)). 

n — 1 

Let j be a temporary sleeper for {x'-''} and let i be a temporary sleeper for {j/*^''}. 
The existence of temporary sleepers was shown in Proposition 14.41 Thus the to- 
tal sum of all logarithmic coordinate losses of y*-'^ at each iteration is at least 
Py^{A'^^\ . . . , A'^'^^), while the ith coordinate of y'-'-' is a sleeper, and the total sum of 
all logarithmic coordinate losses of x^^^^ is at least ;jztPs(^ ; ■ ■ • : ^^'^') while the 
jth coordinate of x^'^ is a sleeper. This will stop the alternating method. Indeed, 
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we repeatedly apply Pk ■ ■ ■ Pi and stop when all coordinates of y^'^ decrease with 
respect to that of y^^\ As y^''\ for / > 1, is a max combination of the columns of 
^(fe)^ by Proposition [4T0l we have that logy^^^^ -logyf^ < ||y(^)|| < \\A''\\ for aU 
t e [n]. Lower bound (|4.16p for the total sum of logarithmic coordinate losses of 
2/(') at each iteration implies that after at most 2{n - l)\\A^'''>\\/p^{A^'^\ . . . , A'^'"')) 
iterations there will be t such that \ogyl — \o%y\ > \\A^'''>\\, if the method does 
not stop, and this contradicts Proposition 14. 101 Hence, after at most that number 
of iterations all coordinates will have to fall in value with respect to the coordi- 
nates of the initial vector. Now, as x^''\ for / > 1, is a min combination of the 
columns of ^C^), by Proposition 14.101 we have that loga;^ — logx^- < ||x'^'|| < 
|^(fe)r|| for all t e [nik] (note that \\A\\ = ||A'^|| for any positive matrix A). 
Using (|4.17p instead of (|4.16p and arguing as above, we obtain the upper bound 
2{mk - l)\\A^'''^'^\\/{-^p^{A''^\. . .,A'-'''>)) on the number of iterations, and this 
proves the claim. D 

If there is more than one matrix with all entries positive, then bound (|4.15p 
can be improved. 

Theorem 4.14. Suppose that A^^^ e R"'""\...,AW g M"''™^. that 
A'-'^'-\ . . . ,74^'"='' have all entries positive, and that span(A(-'-)) n . . . n span(v4*^'^)) — 
{0}. Then after no more than 

(4.18) 2{n - l)mlnmin(p('^-)||, (m,. - l)p(''-)^||)/ps(^^'\ • ■ • , A^''^) 

2—1 

iterations the alternating method will terminate. 



Proof. Applying the argument of Theorem 14.131 and using the fact that ps, 
like ph, is stable under the cyclic permutations of its arguments, we obtain that for 
any t — 1, . . . , s, after at most 

(4.19) l = 2{n- l)min(p(^')||, (m,, - l)p(-*)^||)/ps(AW, . . . , A^)) 

iterations all coordinates of y'-'' ''* have to fall with respect to the coordinates of 
y(i) rt ^ This means that there is a /i < 1 such that y*^'^ ''' < ny^^^ '"* . As all projectors 
are homogeneous and order preserving, we also have that y^' < fJ.y^^' ■ Therefore 
all the coordinates of y^'^^ decrease with respect to that of y'^^\ and hence to that 
of y^^' , and the alternating method stops with negative answer. So the number 
of iterations does not exceed (|4.19p for each rt, and hence it does not exceed the 
minimum of these, which is (|4.18p . D 

Now we show that the techniques developed above apply to the case of integer 
matrices over the max-plus semiring IRmax,+ = (M U {— oo},© — max, (g) = +) 
investigated by Cuninghame-Green and Butkovic 20 . In what follows, we switch 
to the matrix algebra over the max-plus semiring and to the alternating method 
formulated over that semiring. 

First note that if y £ M" is a max-plus or min-plus combination of columns of 
a matrix A E M"x™ with real entries, then ||y|| < P||, where like in ()4.14p but 
without logarithm, the norms are defined by 

(4.20) l|y|| = max(yj - yj), \\A\\ = max{aik - a^k)- 
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Theorem 4.15. Suppose that ^(i) G E"^™i, . . . , A'-'"'^ e M"^""" /lave all entries 
integer. Then after no more than 

(4.21) 2m'inmin((n-l)^— ip('')||,(TO, -1)P«'^||) 

i=l k 

iterations the alternating method will terminate. 

Proof. We are in almost the same situation as in Theorem 14. 141 for all x'-'-*'* 
and y^'^'* there exist temporary sleepers, the norms ||y'-''''*|| do not exceed P'-'*-'!! 
and the norms ||x'^'^''|| do not exceed ||74*^^)-^||. It remains to give bounds for the 
total sum of coordinate losses for cc'''^ and j/'-'-'* at each iteration. As everything 
is integer, the total sum of losses for both x''^'' and y^'^* is not less than 1. The 
multiple ^^ at ||yl'^*)||, which may be important only if k is small, is due to the 
observation that if we apply Pi, . . . ,Pk-i to y^'^ G A^*^^ and do not see any fall 
in coordinates, then y*-'-* is in the intersection and the method immediately stops, 
hence during the run of the algorithm, after at most fc — 1 actions (not k but fc — 1) 
of the sole projectors at least one coordinate of y has to fall. The claim now follows 
by the same argument as in Theorems 14.131 and 14.141 D 



The bounds on number of iterations in [20j . obtained in the case k = 2, are 
in the same vein as (|4.2ip . The only bound on number of iterations in [20j which 
does not depend on the choice of initial vector would read in our terms essentially 
as 2minf^^((m; — 1) max^jda^ |)), where | • | denotes the modulus of an entry. The 
bound of (j4.2ip is expressed in terms of projective norms of rows and columns of 
the matrices, which makes it more precise. 
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